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Abstract 

We study the evolution of quantum fluctuations of a scalar field which 
is coupled to the geometry, in an exponentially expanding universe. We 
derive an expression for the spectrum of intrinsic perturbations, and it 
is shown that the intrinsic degrees of freedom are well behaved in the 
infra-red part of the spectrum. We conclude that quantum fluctuations 
do not change the dynamics of the spacetime in a way which makes its 
evolution non-perturbative and stochastic. This result contradicts previ¬ 
ous derivations which are based on the study of a quantum field on a fixed 
geometry. 

1 Introduction 

It is by now widely believed that the universe has gone through an early pe¬ 
riod of nearly exponential expansion (see e.g. [Q for an overview). Apart from 
solving the horizon, flatness and monopole problems, the assumption of nearly 
exponential expansion has the important implication of producing a nearly scale 
independent spectrum of perturbations, which is in agreement with observations 
Quantum fluctuations of the field which drives inflation, the inflaton, are 
generally believed to be the source of the perturbations which gave rise to the 
observed structure in the universe at late times. An idea which goes further has 
been proposed by Linde ||^ - ||], who argues that quantum fluctuations of the 
inflaton grow in amplitude and decohere during inflation. Due to decoherence, 
the state of the inflaton can be described by a probability distribution, and due 
to quantum fluctuations, the time evolution of the inflaton becomes stochastic. 
Although the mean field of the inflaton tends to evolve to a minimum of its 
potential, the spread of the probability distribution in time allows for a mo¬ 
tion away from the minimum of the potential, with some finite probability. If 
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there is a sufficiently large probability of finding the inflaton at a larger value 
of the potential, then inflation could proceed indefinitely in certain regions of 
the spacetime, which then dominate in terms of the spatial volume due to the 
expansion, while inflation may have come to an end in other regions. The large- 
scale structure of the spacetime can then be expected to be very irregular, hence 
the name ‘stochastic inflation’. In this paper we will reconsider Linde’s idea of 
stochastic inflation, by studying the evolution of a quantized scalar field which 
drives inflation on a perturbed rather than a fixed geometry. The important 
difference between our derivation and the derivations in §- § arises from the 
definition and the interpretation of the perturbations. Since in previous treat¬ 
ments one considered inflaton perturbations on a fixed geometry, it was possible 
to define a probability for finding the inflaton field in a specific state at some 
time. This picture breaks down when one considers quantum fluctuations of the 
inflaton field, which are coupled to the geometry, since then there is no notion 
of a point in space and time at which one can compare perturbations of the 
inflaton field. Instead, one has to deal with the intrinsic perturbations of the 
inflaton and the geometry, on a background geometry which can be taken to 
be a Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. If the intrinsic 
perturbations of the inflaton and the geometry tend to become nonperturba- 
tively large at late times or large distances, then this shows that the spacetime 
evolves away from FLRW, as is the case in the stochastic inflation picture. We 
will show that this does not happen when the coupling between the inflaton and 
the geometry perturbations is taken into account. 

The paper is structured as follows. In section || we discuss the generation of 
long-range quantum correlations of a scalar field, in an unperturbed exponen¬ 
tially expanding FLRW spacetime. It appears that these quantum correlations 
diverge in the limit where the mass of the field tends to zero. This effect can be 
explained from two different viewpoints; in comoving coordinates we find that 
individual quantum fluctuations ‘freeze out’, while in a static coordinatization 
the quantum field appears to ‘heat up’. 

In section js] we discuss the issue of decoherence, and the generation of clas¬ 
sical perturbations from quantum perturbations. 

The quantization of the inflaton and geometry perturbations is considered in 
section^, where we apply Bardeen’s approach [|| to describe the intrinsic degrees 
of freedom of the system. The spectrum of intrinsic quantum perturbations is 
obtained in the idealized case of vanishing inflaton mass and exact exponential 
expansion, and we study the nature of the geometry perturbations at late times 
and large length scales. 

It is shown that the intrinsic deviation of the perturbed spacetime from 
FLRW does not grow at large length scales or at late times. This observation, 
and the observational bound on the amplitude of the perturbations, justify 
our perturbative approach to the problem. We conclude in section || with the 
observation that no nonperturbative effects occur of the kind which could drive 
stochastic inflation. 
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From now on we choose our units so that Ti = c = 8nG = 1, and we let 
Greek indices run from 0 to 3, while Latin indices run from 1 to 3. 


2 The generation of fluctuations 

In this section we will consider the generation of quantum fluctuations in an 
exponentially expanding FLRW spacetime, and its analytic extension which is 
de Sitter spacetime. The use of exactly exponentially expanding FLRW has the 
advantage that there are exact results available (see e.g. |P |^), which is not 
the case for more general inflating spacetimes. Although de Sitter spacetime and 
a nearly exponentially expanding FLRW spacetime are topologically very dif¬ 
ferent, there are good reasons to believe that locally defined physical quantities 
are approximately the same in these two spacetimes. This observation is known 
as the ‘no hair conjecture’ for exponentially expanding spacetimes ||l^. The 
metric of spatially flat exponentially expanding FLRW in terms of comoving 
coordinates reads, 


ds^ = {Hrj) ^ (—d? 7 ^ + JijdxMx^), (1) 

where 77 G R~ denotes the conformal time parameter, and H is the Hubble 
constant. Let us consider a free Klein-Gordon field (/) with a mass m in the ge¬ 
ometry (|^. In the case where the quantum state of the 0-field is invariant under 
spatial translations, the field operator which describes the quantized (()-field can 
be expanded in terms of a basis of modes which are the product of normalized 
plane waves Q, which we define in the appendix, and a time dependent factor 
(j)o which solves the equation 

• • • 7T7 

- 277-Vo + (k^ + (tt= 0 ’ ( 2 ) 

Hr] 

where a dot denotes differentiation with respect to the conformal time 77 , and k 
denotes the comoving wavenumber. In the case where the mass m of the field 
is zero, the solutions of equation (H) are given by, 

00 = c77^iLV(-k77), (3) 
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where denotes the Hankel function of the first kind, and c is a normaliza¬ 
tion constant |Q. At late times, i.e. when 77 I 0, the massless solutions (||) 
approach a nonzero constant, while massive solutions of equation go to zero 

as (— 77 ) 3 ^. Note that although the amplitude of the modes remains large at 
late times when m H, the probability density jo = -^(0*0 — 0*'/') which is 
associated with the modes decreases as fast as . The dynamical picture of the 
quantum fluctuations in an exponentially expanding FLRW spacetime is that 
their amplitude remains large at late times, in spite of exponentially vanishing 
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probability densities. In terms of the quantization which is natural in comoving 
coordinates, there is no change in the occupation number of different modes, 


and for this reason there is no creation of particles |12|. 

Let G'o(x) denote the two point correlation function of the massless (/>-field 
which is assumed to be in the vacuum state which is invariant under spatial 
translations. The correlation function G'o(x) is given by a sum over modes, 
which factorize in terms of the time dependent solutions (/)o, and the spatially 
dependent harmonics Q (see e.g. 


Go(x) = lyQ*(k,x)Q(k,0). 


(4) 


Let Go(k) denote the spatial Fourier transform of Go(x), defined by expression 
in the appendix. From expression it follows directly that 

Go(k) = i|0o(?7,k)p. (5) 

Using expression (^) for (f>Q, and the asymptotic behavior of the Hankel function, 
Hq^\z) ~ z~‘^ for z J, 0, we find the asymptotic expression for Go(k), 

Go(k) ~ k-3, (6) 

for k J, 0. Expression (^) cannot be integrated about k = 0, and therefore 
the correlation function Go(x), which is formally defined by expression (||), 
diverges for arbitrary x. This is the well known infra-red divergence of the 
correlation function for a massless scalar field in an exponentially expanding 
FLRW spacetime |§] - [p]| . 

If one is interested in the time dependence of quantum correlations, then 
equation (^) is somewhat difficult to interpret, since the comoving wavenumber 
k, and the physical wavenumber k, are related by the time dependent expression 
k = Hrjk. In order to describe the time dependence of quantum correlations, it 
is natural to consider the Fourier transform of the correlation function (^) with 
respect to harmonics which are normalized in physical fc-space. From expression 
(H) and the relation ( ^l|) in the appendix, it follows that the Fourier transform 
of Go(x) in physical fc-space is given by 

Go(fc) = {Hp)-^Go{k). (7) 

The asymptotic behavior of Go{k) follows from expression (||) and (@), i.e., 

Go(fc) ^ fc-3, (8) 


for k I 0 and rj constant, and there is no time dependence in this equation. 

The infra-red divergence in equation (^) is related to the global properties 
of the exponentially expanding FLRW geometry. If one assumes that the ex¬ 
ponential expansion has started at a finite time in the past, then one can show 


4 



that Go{k) is well behaved about fc = 0 ||]. This follows e.g. by calculating the 
infra-red contribution to the expectation value of 

{\4>\^)= td^kCik), (9) 

Jo 

where A is an upper boundary on the /c-space integral. It can be shown that 
(|(^P) grows linearly in time when the exponential expansion has started at a 
finite time in the past, and approaches a constant for a positive mass m 

. . . 

A different perspective from which one can look at the quantum correlations 
in an exponentially expanding FLRW spacetime, is the presence of a thermal 
distribution of particles which is natural with respect to a static coordinate sys¬ 
tem. On the analytic extension of an exponentially expanding FLRW spacetime, 
which is de Sitter spacetime, one can choose coordinates so that the metric on 
a section of the spacetime takes the static form 

ds^ = H-^{- sin^ 0df + d0^ + cos^ 6»dO^), (10) 


where 0 G [0,7r/2] is related to the commonly used radial coordinate r by r = 
H~^ COS0, and dfl^ is the surface element on the unit two-sphere. We may now 
define quanta of the 0-field, for which the wave function (j)ujim factorizes in terms 
of a time dependent plane wave ( 2 a;)“ 2 e“*‘^*, a two-sphere harmonic Yim, and 
a function of 0 and Gibbons and Hawking [T^ showed that when the 

quantum state of the field 0 is analytic at the horizon 0 = 0, then quanta with 
wave functions (puiim are present with nonzero occupation numbers. 






( 11 ) 


In a spatially flat FLRW spacetime, a thermal distribution of the form ( 0 ) 
does not give rise to divergent long range correlations, even for m = 0, since 
in this case there is a phase space element of the form d^fc which cancels the 
divergence in expression (H) at w = 0. In the quantization which is based on a 
static coordinatization of de Sitter spacetime, the situation is however different 
since a constant time section in the metric (|l^ ) is one half of a three-sphere. 
The phase space in which the distribution function (0 is evaluated consists 
therefore of the continuous variable uj, and the discrete labels l,m. 

Let us now consider the massless case where m = 0, and u € R'*'. The infra¬ 
red contribution to the vacuum expectation value of |0p, is obtained by taking 
the integral over oj about w « 0, of the occupation number, multiplied by the 
squared modulus of the wave function. Indeed, the occupation number rii^oo and 
the square of </>ujoo diverge as near tu = 0, and their formal integral over aj 
diverges. Recall that in the quantization which is based on comoving coordinates 
we found an infra-red divergent vacuum expectation value of |0p, which is due 
to the ‘freezing’ of zero point fluctuations, while no particles are created. In the 
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quantization which is based on static coordinatization, the vacuum expectation 
value of \(j>\’^ is also infra-red divergent, but this time it is due to the ’heating’ 
of the quantum field on spatially bounded sections. 

3 Classical and quantum perturbations 

In the previous section, we discussed the generation of quantum correlations 
in an exponentially expanding FLRW spacetime. In the quantization which 
is based on comoving coordinates, we found quantum fluctuations with large 
amplitudes, although in this description the occupation number remains zero for 
each mode [|l^ . The assumption that large amplitudes imply large occupation 
numbers seems to play a role in a number of references §- § , where it is used 
to argue that long range correlations become classical in some sense. 

The question how quantum correlations at early times evolve into classical 
perturbations at late times, appears to receive little attention in the literature, 
although efforts towards a more satisfactory description have been made (see 
e.g. Q - Q). Without going into much detail, let us go through some of the 
ideas which play a role in the transition from quantum to classical perturba¬ 
tions. What appears to be crucial, is the decoherence of inequivalent histories 
[ p0| . In the context of the discussion which is given in the previous section, 
a history 4>{x) is a configuration of the (/)-field which is given over the entire 
spacetime. In physically relevant situations, it is often the case that only some 
degrees of freedom of the (/)-field are accessible to observation, and it is natural 
to consider the equivalence classes of histories which cannot be distinguished 
by observation. While the state of the 0-field is generally given by a density 
matrix, a ‘reduced density matrix’ can be obtained by projecting the density 
matrix onto sets of inequivalent histories. Decoherence occurs when the reduced 
density matrix is approximately diagonal, and in this case the diagonal elements 
can be interpreted as approximate probabilities for inequivalent histories. When 
two or more sets of histories are mutually exclusive at a classical level (e.g., as 
in Schrodinger’s experiment where there are two states of a cat in a box), then 
they must have decohered to a high degree of accuracy, in order to account for 
the absence of observed interference effects. 

Note that there is a considerable amount of freedom to choose a criterion by 
which histories are considered to be equivalent. Which criterion is natural in a 
certain situation, depends on the type of measurement which is conducted. 

As an example, let us consider again the case of a free scalar field 0, which is 
in the translation invariant vacuum state in an exponentially expanding FLRW 
spacetime. A geodesic observer in this spacetime perceives an event horizon at 
the Hubble distance , and can only respond to quanta which originate from 
somewhere within this event horizon. One can show that a geodesic observer 
responds to quanta which have a positive frequency with respect to the observers 
proper time (see e.g. [^). These quanta are naturally defined in terms of a 
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quantization which is based on a static coordinate system of the form 
where the coordinates are chosen so that the center of the coordinate system at 
6 * = f coincides with the worldline of the geodesic observer. When expressed 
in terms of the quantization which is based on static coordinates, the state 
of the 0-field is described by a reduced density matrix, which appears to be 
completely diagonal [l^ . Decoherence is in this case obtained by summing 
over the degrees of freedom of the 0-field which are outside the event horizon. 
However, the notion of an event horizon is observer dependent, and the reduced 
density matrix is in this case only natural to describe the measurement of one 
specific geodesic observer. 

In the case of the early universe, there are no observers which can absorb 
quanta of a scalar field, and one has to consider what one could call a ‘measure¬ 
ment situation’ [Q . Essentially this means that one introduces a coupling 
of the 0-field with itself or other fields (for instance by adding a cubic interaction 
term in the Lagrangian, or by considering the nonlinear coupling between the 
0-field and gravitational degrees of freedom). The effect of an interaction term 
is that quantum fluctuations at large length scales (compared to the Hubble 
radius H~^), are now coupled to quantum fluctuations at small length scales. 
By summing over perturbations of the 0-field with a wavelength which is small 
compared to the Hubble radius, one obtains a reduced density matrix which 
describes the decohered large scale fluctuations of the 0-field [p^-[l^. 

Let us now assume that decoherence occurs, so that one can define a proba¬ 
bility distribution P(0, t) which describes the probability P of finding the quan¬ 
tum field with a value 0 at a point in the spacetime with time t. In Linde’s 
approach to stochastic inflation §-i, the time evolution of the probability 
distribution P(0, t) is described by a diffusion equation, 

dP{(t>,t) 

dt 02^ ’ ^ ^ 

where I? is a positive constant. The spread of the quantum field 0 about its 
mean value (0), is given by, 

(i02|)=y'd0i0-(0)pp(0,t). (13) 

Recall that in section |2| we calculated (|0^|) for the case of a free field in expo¬ 
nentially expanding FLRW. Although the expression was found to be divergent 
in the case of exact exponential expansion, a linear growth in time can be found 
when the exponential expansion has started at a finite time in the past |^. 
This result can be used to determine the value of the diffusion constant D in 
expression (0) (see e.g. [Q). 

An important limitation of the derivation which we sketched above, is that 
the probability distribution P{(j),t) is defined for an FLRW geometry. However, 
since we are interested in perturbations at length scales large compared to the 
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Hubble radius, there is no good reason why one can ignore the coupling between 
the inflaton perturbations and the geometry 0. But if one does consider the 
coupling between the inflaton perturbations and the geometry, then we do not 
have a single geometry. Instead, for every configuration of the inflaton, one has 
a different geometry, and there is no notion of a single point in space and time 
for which one can define a probability t) for finding the field with a value 
at time t. Indeed, one can show that quantum fluctuations of the inflaton field 
can be interpreted equally well as quantum fluctuations of only the metric, in a 
coordinate system where there is no inflaton perturbation at all. 

In the following section, we derive the spectrum of intrinsic perturbations 
of the inflaton and the geometry which are generated during inflation. 


4 Quantum fluctuations and stochastic inflation 

In the previous sections, we studied the amplification of quantum fluctuations of 
a scalar field in an exponentially expanding FLRW spacetime. The correlation 
function was found to be infra-red divergent in the case of a massless scalar 
field, while a finite result was found for a positive scalar field mass or when 
we assumed that the exponential expansion had started at a finite time in the 
past. In the physically more realistic case of a universe where approximate ex¬ 
ponential expansion is driven by a scalar field with a potential term, the scalar 
field mass will generally be nonzero and time dependent, and the expansion will 
tend to slow down gradually before inflation ends. However, since our aim is 
to investigate the possibility of an unbounded growth of intrinsic perturbations 
during inflation, we will neglect the mass of the inflaton and the slowing down of 
the expansion while the inflaton moves down the potential. In our perturbative 
approach to describe the evolution of quantum fluctuations of the inflaton and 
the geometry during inflation, we will therefore adopt a background geometry 
of the form ©■ As is usual in cosmology, we define a perturbation of a physical 
quantity as the difference of the same physical quantity, evaluated at corre¬ 
sponding points in a perturbed and a background spacetime. What one calls 
a perturbation therefore depends on a mapping between points in a perturbed 
and a background spacetime, which is called a gauge. In the following we will 
use Bardeen’s formalism [Q to deal with the degrees of freedom which, to linear 
order, do not depend on the choice of gauge. In our description of a scalar field 
in a perturbed spatially flat FLRW spacetime, it will be sufficient to consider 
only scalar perturbations, and the metric can be expanded as 

[-^°<5°(1 + 2A)Q + {Sl6l + 6^X)BQ,, ++ 2iJiQ) 

k 


-\-2HTh^^hi, (k ^Q-ij + 


(14) 
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where ; i denotes the derivative with respect to the comoving spatial coordinate 
x% hii := is the projection operator onto spatial hypersurfaces of 

constant time in the background, and the harmonics Q = Q(x, k) are solutions 
of the scalar Helmholtz equation which are normalized according to equation 
( p4[ ) in the appendix. We assume that inflation is driven by a scalar field, which 
we call Ip in order to make a distinction with the held p which was dehned on a 
hxed geometry in the previous sections. The Lagrangian density of the held ip 
is given by 

(15) 

where the potential term Vpip) is bounded from below. It will be useful to write 
the held ip as the sum of a perturbation Sippx.,t), and the mean value ipopt), 
which is dehned as the average of ip^ evaluated on a spatial hypersurface of 
constant time in the background. 

The mass of the held ip is dehned as the curvature of the potential, i.e. mP := 
^d'^V(ip)/dip'^, which is evaluated at ip = ipQ. Nearly exponential expansion 
occurs when the kinetic term in the Lagrangian ( p^ ) is small compared to the 
potential term (so that Tjf « ^S^Vpipo)), and V{ipo) is approximately constant 
on a timescale of the expansion time. It can be shown that nearly exponential 
expansion can occur for a large class of potentials, including potentials of the 
polynomial form V{ip) = with A G i?+ and n G Z+ (see e.g. |^). 

In the following, we will denote perturbations by their Fourier components, 
which we dehned in the appendix, and we drop the argument k, unless there is 
a risk of confusion. 

A gauge invariant variable which is useful in the description of perturbations 
is given by. 

Pm ■■= Sip+ ^ {Hl (16) 

and pm can be interpreted in terms of the three-curvature of the constant-';/’ 
hypersurfaces, 

'Wonstanti/' “ ^ 
where we assume that pQ ^ Q. 

Copying Bardeen’s notation, we dehne Cm as the fractional energy density 
perturbation in the comoving time-orthogonal gauge, i.e., 

Em ■■= ST°/T°, (18) 

evaluated in the gauge where 6p = B = 0. 

A variable related to Em is Om, which is dehned as the fractional lapse 
function perturbation, evaluated in the comoving time-orthogonal gauge, i.e. 
am := A in the gauge where Sp = B = 0. From the expression for the energy 
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density of the scalar field in the comoving time-orthogonal gauge, i.e. Tq = 
-I- it follows that 


Cm — (t “f ^)crm; 


(19) 


. 2 

where uj := —l-\-2l{l + V {ij)) /^o )■ The fractional energy density perturbation 
Cm is coupled to the entropy perturbation rje, defined as the difference between 
the fractional isotropic pressure in the perturbed spacetime, and the fractional 
isotropic pressure at a point in the background spacetime with the same en¬ 
ergy density. The entropy perturbation rje, and the fractional lapse function 
perturbation Um, can be shown to be proportional 0, 




UJ 

(l+u;)(c2-l)^' 


( 20 ) 


where we specialized Bardeen’s equation (5.20) to the case with scalar field 
matter, and the speed of sound can be expressed in terms of background 
variables. Hence, Cm, Om and rje all represent the same physical perturbation, 
but unlike Cm and Om, the interpretation of rje is not related to the properties 
of some collection of spatial hypersurfaces in the perturbed spacetime, and it is 
therefore an intuitively clear measure of an intrinsic perturbation. 

The gravitational and infiaton action has been expanded to second order in 
terms of the scaled perturbation variable w := S/pm and background quantities, 
by Deruelle et al. [|2|. Recall that the variable prm defined by expression 
(p^, is only to linear order invariant under a gauge transformation. The square 
of a first order gauge invariant variable, is however gauge invariant to second 
order. An expansion of the action which contains terms of linear and quadratic 
order in is not in general gauge invariant to second order, unless the linear 
terms in the expansion vanish. Note that an expansion of the action in terms 
of the perturbation variable is obtained by integrating the expansion of 
the Lagrangian dil) and the volume element, over the background spacetime. 
Indeed, terms which are linear in pm vanish in this integration, due to the 


orthogonality relation for the harmonics Q (equation (34) in the appendix). 

The field equation for w appears to be of a simple Klein-Gordon type with 
a time dependent mass term. 


- 2?7 ‘^)w(k) = 0 


( 21 ) 


which has the solutions 


a;(k) = (-,yk). 


( 22 ) 


where denotes the Hankel function of the first kind [^. Note that equation 
(13) is derived from an expansion of the action of the infiaton field and the 
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geometry, in terms of the gauge invariant variable w | |^ . The coupling between 
the inflaton and the geometry is therefore accounted for by the requirement that 
the solutions ( |2^ ) extremize the combined gravitational and inflaton action. 

The field operator which is associated with the second quantized variable w 
can be expanded in terms of the solutions (|2^, and the two-point correlation 
function for the field w takes the usual form [|12[, 

Go(k) = i |zc(k)|^ (23) 

where we assumed that the field is in the vacuum state which is natural in 
comoving coordinates. In the following we assume that the quantum state of 
the u>-field has decohered (see section^), in which case expression (|^ ) can be 
interpreted as the sum over decohered configurations {rc} of the square of the 
perturbation component w(k) evaluated for {w}, multiplied by the probability 
of finding the configuration {w}. 

The asymptotic behavior of Go(k) near k = 0 follows from expression (p^, 
( p^ and the expansion of the Hankel function, H^\z) ^ z~'^ for z J, 0. We find 

Go(k) - ry-2k-3, (24) 

for k I 0 and rj constant. Assuming decoherence, expression ( p4|) describes the 
asymptotic behavior of the classical perturbation component ■u;(k) for k J, 0, 

w(k) ~ 7;-ik-5, (25) 

where here and in the following we neglect a probabilistic factor of order one on 
the right-hand side of the ^ symbol. 

4.1 Energy and lapse function perturbations 

Let us now determine the spectrum of fractional energy density and lapse func¬ 
tion perturbations. The fractional energy density perturbation Cm has been 
evaluated in terms of w in Q , 

kH 

= 2—[w - ilS'w], (26) 

where k is a constant which depends on the background geometry. Using ex¬ 
pression (p^ ) for w, and the differentiation property of the Hankel function, 
{d/<lz)Hq{z) + {q/z)Hq{z) = ilg_i(z), we find, 

^ ky=i) (27) 

where we used that the scale factor S equals {Hr])~^ in the case of exponential 
expansion. Using the asymptotic expression for the Hankel function, Hq^\z) ^ 
for z i 0, yields. 

Cm ~ CXm ^ Tj k^ . (28) 
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The time dependence in equation ( p^ is a consequence of the Fourier decom¬ 
position of the perturbations in comoving rather than physical fc-space. Let Zm 
and (Xm denote the physical fc-space Fourier components of the fractional energy 
density and lapse function perturbations, which we define by expression ( |40| ) in 
the appendix. According to expression ( pl| ) in the appendix, the perturbation 
components in physical and comoving fc-space are related by and 

dm = S^Um- Using this relation, and expression (p8|), yields the asymptotic 
expression, 

(29) 

where k = S'“^k denotes the physical wavenumber, and there is no time depen¬ 
dence at the right-hand side of equation (|^). Expression ( p^ ) shows that the 
infra-red behavior of the fractional density perturbation and lapse function is 
regular, and integrable about k = 0. 

Note that expression (|^) determines the asymptotic behavior of the spec¬ 
trum of fractional energy density perturbations, but not its magnitude. Obser¬ 
vations of perturbations at the time of last scattering |^, put an upper limit 
on Cm of the order of 10“® for wavenumbers small compared to the value of the 
Hubble parameter at that time. Since Cm grows approximately as the square of 
the scale factor for the growing mode and small wavenumbers |Q , the fractional 
energy density perturbations must have been many orders of magnitude smaller 
than 10“® at the time when inflation came to an end. 

4.2 Curvature perturbations 

In this subsection, we investigate the spectrum of curvature perturbations of the 
hypersurfaces on which the inflaton held tp is constant. According to expression 
(0) in the appendix, the perturbation components in physical and comoving k- 
space are related by (pm = Sicpm- Using expression (|^), we hnd the asymptotic 
behavior of (pm for fc J, 0, 

(pm = S^(pm k~i, (30) 

and there is no time dependence in equation (^0[). 

From expressions (0 and ( |25| ) we derive the expression for the Fourier 
component of the spatial curvature perturbation, with respect to the physical 
wavenumber fc, 

"^constantj/> 

for fc I 0. 

We will now address the question whether the geometry of the constant-')/' 
hypersurfaces develops increasingly strong inhomogeneities at large length scales 
or at late times. More generally, one would like to know whether our perturbed 
FLRW spacetime remains everywhere close to FLRW in some well dehned sense. 

The question whether an inhomogeneous spacetime is close to FLRW in some 
objective sense, is known as the averaging problem in cosmology. This problem 
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appears to be surprisingly nontrivial, even in the linearized case, since it involves 
comparing tensors in different spacetimes, in a way which is independent of the 
choice of coordinates, gauge, and the averaging prescription (see e.g. and 
references therein). 

In the following we will sidestep the difficulties involved with a precise for¬ 
mulation of the averaging problem, by considering the scaling behavior of the 
spectrum of perturbations. Without presenting an objective criterion for a 
spacetime being close to an FLRW spacetime, we determine whether intrinsic 
perturbations grow or decrease under a change of the length scale. This allows 
us to compare the warping of our universe at superhorizon scales and at length 
scales where observations can be made, e.g., on the intersection of our past light 
cone with the surface of last scattering. If perturbations tend to grow at large 
length scales, then our perturbative treatment can be expected to break down 
at superhorizon scales, which is in agreement with the picture of a stochastically 
inflating universe which is strongly warped at superhorizon scales [^. 

On the contrary, a decreasing or constant spectrum of perturbations implies 
that our perturbative approach holds on superhorizon scales, given that it holds 
at subhorizon scales. This will then give us an improved estimate of the effect 
of quantum fluctuations on the global structure of the universe. 

We define a new length scale I by choosing a new unit of physical length 
which is a G R times the former unit of length. In terms of the new length 
scale, a physical quantity q with dimension transforms to g = a~^q. As we 
derived in the appendix, the Fourier transform q{k) of a quantity q(x) transforms 
under scale transformations as a half-density of this quantity in fc-space, i.e. 
q{k) = q{k), where k = ak. 

The scaled spatial curvature perturbation TZ^k) therefore takes the form 

^constantcji ~ ^constant</i> (^2) 

where we used that the spatial curvature has the dimension 

Using expression ( |3l|) for the spectrum of three-curvature perturbations, and 
the scale transformation property (|^), we find, 

^constant<^'^ ® ^ ^ ^ ( 33 ) 

Comparing equations (^) and (^), we find that the spectrum of curvature 
perturbations of the constant-'i/f hypersurfaces remains constant under a change 
of the length scale. 

Summarizing the results derived in this subsection, we found that perturba¬ 
tions in the spatial curvature are regular about fc = 0, and the perturbations do 
not grow at late times or at large length scales. 

4.3 Interpretation 

Let us now address the question why the gauge invariant perturbations TZ and 
dm do not show the same infra-red divergence as we found for a quantum field (j) 
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on a fixed geometry. At first sight, this might be surprising, since expressions (||) 
and ( p^ ) show that the gauge invariant perturbations (pm , and the perturbations 
0 on a fixed geometry, diverge with the same power when fc J, 0. The essential 
difference between the perturbations <pm and (p, is their interpretation in terms of 
physical quantities. Note that a perturbation of the field on a fixed geometry is 
proportional to a perturbation in the energy density, which follows by expanding 
the energy density Tq — ^{d^(pdop+V{p)) about the background value p = 
where we use that V'{p)\^Q 0 as long as the scalar field has not reached 

a minimum of the potential. On the contrary, it follows from equation 
and (30) that the gauge invariant amplitude pm acts as a potential for the 
fractional lapse function and energy density perturbations in the limit when 
fc I 0. Similarly, expression (0) shows that pm acts as a potential for the 
spatial curvature perturbation TZ. The perturbation pm may therefore diverge 
as for fc I 0, while the spatial curvature TZ, and the fractional lapse function 
and energy density perturbations am and are well behaved when fc | 0. 

One may question whether the time evolution of the perturbations in the 
spatial curvature and the fractional lapse function can be described by a diffusion 
equation of the form di)- For a massless quantum field 0 on a fixed geometry, 
this method is natural since the expectation value of \p\'^, which according to 
equation ( ]l^ equals to the standard deviation of the probability distribution 
P{p,t), appeared to grow linearly in time. In this section we showed that 
the spectrum of fractional lapse function and spatial curvature perturbations 
is integrable about fc = 0, and constant in time. The probability distribution 
which describes the decohered fractional lapse function and spatial curvature 
perturbations, is therefore time independent, and it does not evolve according 
to the diffusion equation 


5 Conclusions 

In the previous section we derived expressions for the spectrum of the fractional 
lapse function and spatial curvature perturbations which are generated during 
exponential expansion. We found that intrinsic perturbations of the scalar field 
and the geometry do not grow at late times or at large length scales, during 
exponential expansion. The constancy of the spectrum of perturbations in space 
and time, and the observational bounds on the magnitude of the perturbations, 
justify our perturbative approach. Our result contradicts the assumption that 
quantum fluctuations grow nonperturbatively during inflation, which underlies 
the idea of stochastic inflation. This indicates that an inflating universe is not 
likely to develop a highly irregular structure at superhorizon scales, unless these 
irregularities are present in the initial conditions. 
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7 Appendix 

In this appendix, we discuss the scaling behavior of the flat space Fourier trans¬ 
form. Let us define the solutions Q(k, x), which satisfy the scalar Helmholtz 
equation, i.e. Qll = jkpQ, where ; i denotes the covariant derivative with respect 
to X*, and the solutions Q are normalized according to 


J d^xQ*(k, x)(5(k', x) =(5^(k—k'). 


(34) 


The Helmholtz equation and the normalization condition (34) are satisfied by 
distributions of the form (3(k, x) = where the normalization factor N is 

defined as the distribution which is constant for all x, and which is normalized 
as. 


d^xA^^ = { 2 tt ) 


-3 


(35) 


The Fourier transform of a function /(x) with respect to the comoving wavenum¬ 
ber k is given by, 

/(k) = y d^xQ(k,x)/(x). (36) 

Let us now consider how /(k) transforms under scale transformations. We 
implement a scale transformation by defining a new unit of length £, which 
equals a G R times the former unit of length. In terms of the new length scale, 
we define solutions Q of the Helmholtz equation, which are normalized by the 
condition 

J d^xQ*{k,x)Q{k\x) =S^{k-k'). (37) 

A basis of solutions Q which satisfy the normalization condition (|^) is given by 
Q{k,x) = Ne''^^, where N is defined as the constant distribution which satisfies 
the normalization condition 


J d^xN^ = (27r)-3. (38) 

It follows from the definitions ( ^^ and (|^) that the normalization constants N 
and N are related by N = aiN, which implies 

Q(k,x) = a'^Q{k,x), (39) 


15 



where x = a and k = ak. The Fourier transform of a function with respect 
to the barred wavenumber k is defined by 

f{k) := J d^x Q(k,x)f{x), (40) 

and from the relation (|^) and the definitions (^) and (^) it follows that 

f{k)=a-if{k), (41) 

where k = ak. Expression (iiD shows that the Fourier transform of a function 
f{x) transforms under a scale transformation as a half-density in A:-space. In 
some calculations it will be useful to express the Fourier transform with respect 
to the comoving wavenumber (]3^), in terms of the Fourier transform with respect 
to the physical wavenumber. Since the comoving wavenumber and the physical 
wavenumber are related by a time dependent scale transformation k = S^^k, 
this is simply a special case of the scale transformation discussed above, and 
f{k) and /(k) are related by /(fc) = Sif{k), where k = ^“^k. 
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